Bioinformatics for RNA‐Seq Data Analysis

While RNA sequencing (RNA‐seq) has become increasingly popular for transcriptome profiling, the analysis of the massive amount of data generated by large‐scale RNA‐seq still remains a challenge. RNA‐seq data analyses typically consist of (1) accurate mapping of millions of short sequencing reads to a reference genome, including the identification of splicing events; (2) quantifying expression levels of genes, transcripts, and exons; (3) differential analysis of gene expression among different biological conditions; and (4) biological interpretation of differentially expressed genes. Despite the fact that multiple algorithms pertinent to basic analyses have been developed, there are still a variety of unresolved questions. In this chapter, we review the main tools and algorithms currently available for RNA‐seq data analyses, and our goal is to help RNA‐seq data analysts to make an informed choice of tools in practical RNA‐seq data analysis. In the meantime, RNA‐seq is evolving rapidly, and newer sequencing technologies are briefly introduced, including stranded RNA‐seq, targeted RNA‐seq, and single‐cell RNA‐seq

2008 Inter-laboratory Comparison Study of a Reference Material for Nutrients in Seawater

Autoclaved natural seawater collected in the North Pacific Ocean was used as a reference material for nutrients in seawater (RMNS) during an inter-laboratory comparison (I/C) study conducted in 2008. This study was a follow-up to previous studies conducted in 2003 and 2006. A set of six samples was distributed to each of 58 laboratories in 15 countries around the globe, and results were returned by 54 of those laboratories (15 countries). The homogeneities of samples used in the 2008 I/C study, based on analyses for three determinants, were improved compared to those of samples used in the 2003 and 2006 I/C studies. Results of these I/C studies indicate that most of the participating laboratories have an analytical technique for nutrients that is sufficient to provide data of high comparability. The differences between reported concentrations from the same laboratories in the 2006 and 2008 I/C studies for the same batch of RMNS indicate that most of the laboratories have been maintaining internal comparability for two years. Thus, with the current high level of performance in the participating laboratories, the use of a common reference material and the adaptation of an internationally accepted nutrient scale system would increase comparability among laboratories worldwide, and the use of a certified reference material would establish traceability. In the 2008 I/C study we observed a problem of non-linearity of the instruments of the participating laboratories similar to that observed among the laboratories in the 2006 I/C study. This problem of non-linearity should be investigated and discussed to improve comparability for the full range of nutrient concentrations. For silicate comparability in particular, we see relatively larger consensus standard deviations than those for nitrate and phosphate

A Timescale for Evolution, Population Expansion, and Spatial Spread of an Emerging Clone of Methicillin-Resistant Staphylococcus aureus

Due to the lack of fossil evidence, the timescales of bacterial evolution are largely unknown. The speed with which genetic change accumulates in populations of pathogenic bacteria, however, is a key parameter that is crucial for understanding the emergence of traits such as increased virulence or antibiotic resistance, together with the forces driving pathogen spread. Methicillin-resistant Staphylococcus aureus (MRSA) is a common cause of hospital-acquired infections. We have investigated an MRSA strain (ST225) that is highly prevalent in hospitals in Central Europe. By using mutation discovery at 269 genetic loci (118,804 basepairs) within an international isolate collection, we ascertained extremely low diversity among European ST225 isolates, indicating that a recent population bottleneck had preceded the expansion of this clone. In contrast, US isolates were more divergent, suggesting they represent the ancestral population. While diversity was low, however, our results demonstrate that the short-term evolutionary rate in this natural population of MRSA resulted in the accumulation of measurable DNA sequence variation within two decades, which we could exploit to reconstruct its recent demographic history and the spatiotemporal dynamics of spread. By applying Bayesian coalescent methods on DNA sequences serially sampled through time, we estimated that ST225 had diverged since approximately 1990 (1987 to 1994), and that expansion of the European clade began in 1995 (1991 to 1999), several years before the new clone was recognized. Demographic analysis based on DNA sequence variation indicated a sharp increase of bacterial population size from 2001 to 2004, which is concordant with the reported prevalence of this strain in several European countries. A detailed ancestry-based reconstruction of the spatiotemporal dispersal dynamics suggested a pattern of frequent transmission of the ST225 clone among hospitals within Central Europe. In addition, comparative genomics indicated complex bacteriophage dynamics

Performance of non-invasive tests and histology for the prediction of clinical outcomes in patients with non-alcoholic fatty liver disease: an individual participant data meta-analysis

BackgroundHistologically assessed liver fibrosis stage has prognostic significance in patients with non-alcoholic fatty liver disease (NAFLD) and is accepted as a surrogate endpoint in clinical trials for non-cirrhotic NAFLD. Our aim was to compare the prognostic performance of non-invasive tests with liver histology in patients with NAFLD.MethodsThis was an individual participant data meta-analysis of the prognostic performance of histologically assessed fibrosis stage (F0–4), liver stiffness measured by vibration-controlled transient elastography (LSM-VCTE), fibrosis-4 index (FIB-4), and NAFLD fibrosis score (NFS) in patients with NAFLD. The literature was searched for a previously published systematic review on the diagnostic accuracy of imaging and simple non-invasive tests and updated to Jan 12, 2022 for this study. Studies were identified through PubMed/MEDLINE, EMBASE, and CENTRAL, and authors were contacted for individual participant data, including outcome data, with a minimum of 12 months of follow-up. The primary outcome was a composite endpoint of all-cause mortality, hepatocellular carcinoma, liver transplantation, or cirrhosis complications (ie, ascites, variceal bleeding, hepatic encephalopathy, or progression to a MELD score ≥15). We calculated aggregated survival curves for trichotomised groups and compared them using stratified log-rank tests (histology: F0–2 vs F3 vs F4; LSM: 2·67; NFS: 0·676), calculated areas under the time-dependent receiver operating characteristic curves (tAUC), and performed Cox proportional-hazards regression to adjust for confounding. This study was registered with PROSPERO, CRD42022312226.FindingsOf 65 eligible studies, we included data on 2518 patients with biopsy-proven NAFLD from 25 studies (1126 [44·7%] were female, median age was 54 years [IQR 44–63), and 1161 [46·1%] had type 2 diabetes). After a median follow-up of 57 months [IQR 33–91], the composite endpoint was observed in 145 (5·8%) patients. Stratified log-rank tests showed significant differences between the trichotomised patient groups (p<0·0001 for all comparisons). The tAUC at 5 years were 0·72 (95% CI 0·62–0·81) for histology, 0·76 (0·70–0·83) for LSM-VCTE, 0·74 (0·64–0·82) for FIB-4, and 0·70 (0·63–0·80) for NFS. All index tests were significant predictors of the primary outcome after adjustment for confounders in the Cox regression.InterpretationSimple non-invasive tests performed as well as histologically assessed fibrosis in predicting clinical outcomes in patients with NAFLD and could be considered as alternatives to liver biopsy in some cases

Brown fat determination and development from muscle precursor cells by novel action of bone morphogenetic protein 6.

Brown adipose tissue (BAT) plays a pivotal role in promoting energy expenditure by the virtue of uncoupling protein-1 (UCP-1) that differentiates BAT from its energy storing white adipose tissue (WAT) counterpart. The clinical implication of "classical" BAT (originates from Myf5 positive myoblastic lineage) or the "beige" fat (originates through trans-differentiation of WAT) activation in improving metabolic parameters is now becoming apparent. However, the inducers and endogenous molecular determinants that govern the lineage commitment and differentiation of classical BAT remain obscure. We report here that in the absence of any forced gene expression, stimulation with bone morphogenetic protein 6 (BMP6) induces brown fat differentiation from skeletal muscle precursor cells of murine and human origins. Through a comprehensive transcriptional profiling approach, we have discovered that two days of BMP6 stimulation in C2C12 myoblast cells is sufficient to induce genes characteristic of brown preadipocytes. This developmental switch is modulated in part by newly identified regulators, Optineurin (Optn) and Cyclooxygenase-2 (Cox2). Furthermore, pathway analyses using the Causal Reasoning Engine (CRE) identified additional potential causal drivers of this BMP6 induced commitment switch. Subsequent analyses to decipher key pathway that facilitates terminal differentiation of these BMP6 primed cells identified a key role for Insulin Like Growth Factor-1 Receptor (IGF-1R). Collectively these data highlight a therapeutically innovative role for BMP6 by providing a means to enhance the amount of myogenic lineage derived brown fat

Additional file 1: Table S1. of QuickMIRSeq: a pipeline for quick and accurate quantification of both known miRNAs and isomiRs by jointly processing multiple samples from microRNA sequencing

Human mature miRNAs in miRBase Release 21 with identical sequences. Table S2. Human hairpins in miRBase Release 21 with identical sequences. Table S3. Pairs of miRNAs that are reverse complementary to each other in human miRBase Release 21. Table S4. Top 10 miRNAs with large differences in miRNA quantification between stranded and non-stranded mapping modes. Table S5. Distribution of 5′ and 3′ end offsets of unique miRNA reads in GSE64977. Figure S1. Top panel: All of the miRNAs in the alignment have the same mature sequence (highlighted in gray), but originate from different genes as evidenced by the differences in the pre-miRNA sequences. Bottom panel: miRNA genes found in a cluster on human chromosome 19. Figure S2. Protocol of isomiR quantification. Figure S3. Scatter plots of miRNA quantification results by miRge for samples SRR1759212 SRR1759213, SRR1759214 and SRR1759215. The same dataset were analyzed with and without incorporation of the strand information, respectively. Figure S4. Breakdown of mapped miRNA reads into perfect and mismatch categories. Figure S5. Comprehensive annotation of miRNA-seq reads. The summary plot provides an overview of the distribution of annotated reads in all five annotated RNA categories for each sample. Figure S6. Summary report for adapter trimming. Figure S7. Read length distributions for samples SRR1759212, SRR1759213, SRR1759214, and SRR1759215 in the GSE64977 miRNA-seq dataset. Figure S8. Variations at 5′ and 3′ ends of miRNA reads. Figure S9. The comparison of QuickMIRSeq with miRge. (PDF 1088 kb

Characterizing Early Psychosocial Functioning of Parents of Children with Moderate to Severe Genital Ambiguity due to Disorders of Sex Development

PurposeWe examined the psychosocial characteristics of parents of children with disorders of sex development at early presentation to a disorders of sex development clinic. Parental anxiety, depression, quality of life, illness uncertainty and posttraumatic stress symptoms were assessed. Additionally we evaluated the relationship of assigned child gender to parental outcomes.Materials and methodsA total of 51 parents of children with ambiguous or atypical genitalia were recruited from 7 centers specializing in treatment of disorders of sex development. At initial assessment no child had undergone genitoplasty. Parents completed the Cosmetic Appearance Rating Scale, Beck Anxiety Inventory, Beck Depression Inventory, SF-36, Parent Perception of Uncertainty Scale and Impact of Event Scale-Revised.ResultsA large percentage of parents (54.5%) were dissatisfied with the genital appearance of their child, and a small but significant percentage reported symptoms of anxiety, depression, diminished quality of life, uncertainty and posttraumatic stress. Few gender differences emerged.ConclusionsAlthough many parents function well, a subset experience significant psychological distress around the time of diagnosis of a disorder of sex development in their child. Early screening to assess the need for psychosocial interventions is warranted